Critical Trap Aspect Ratios For Dipolar BEC 
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We show that there exists critical trap aspect ratios for a trapped Bose-Einstein condensate with 
dipole-dipole interactions. We discuss the role of critical trap aspect ratios on both the critical angu- 
lar velocity above which a vortex is energetically favorable and the precession velocity of an off-axis 
vortex in TF regime. We show that the stability diagram for a purely dipolar gas depends crucially 
on the critical trap aspect ratios both for positive and negative strength of dipolar interaction. We 
prove that the critical points can be changed by an anharmonic potential. 
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I. INTRODUCTION 



The ultracold gases with dipole-dipole interaction 
became the subject of intensive theoretical studies 
after the first gaseous condensate was produced. The 
experimental realization of Bose-Einstein condensation 
of ^^Cr atoms in 2005 has raised the interest in study 
of BEC with nonlocal dipole-dipole interactions 
This nonlocal character of the potential has remarkable 
consequences for the physics of rotating dipolar gases in 
TF limit. It is believed that, in axially symmetric traps 
with the axis along the dipole orientation, the critical 
angular velocity, above which a vortex is energetically 
favorable, is decreased due to the dipolar interaction 
in oblate traps, and increased in prolate traps 043] ■ 
It was noted in Q that the effect of the dipole-dipole 
interaction is the lowered (raised) precession velocity of 
an off-center straight vortex line in an oblate (a prolate) 
trap. 

The nonlocal character of dipolar potential is also crucial 
for the stability properties of dipolar gases. It was shown 
that there exists a critical aspect ratio in the stability 
diagram of a purely dipolar condensate @ . 
In this paper, we show that the crossover between a 
reduced and an increased critical angular velocity above 
which a vortex is energetically favorable occurs at a 
critical aspect ratio in TF regime. We discuss the role of 
the critical aspect ratio on the precession velocity of an 
off-center straight vortex in TF limit. In the Gaussian 
limit, we show that there exists two critical aspect ratios 
for a purely dipolar condensate. We investigate their role 
on the stability diagram. We discuss that an additional 
anharmonic potential changes the critical trap aspect 
ratios. 

The paper is organized as follows: The following section 
reviews two main methods to calculate the dipolar 
potential interaction. Section III finds a critical trap 
aspect ratio in TF limit. Section IV investigates critical 
trap aspect ratio in the stability diagram and discusses 
the effect of additional quartic potential on the critical 
aspect ratios. 



II. DIPOLAR INTERACTION 

Consider a BEC of N particles with mass m and mag- 
netic dipole moment fi oriented in the z-direction by a 
sufficiently large external field. At sufficiently low tem- 
peratures, the description of the ground state of the 
condensate is provided by the solution of the Gross- 
Pitaevskii (GP) equation 
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mean field term due to dipole-dipole interactions, Cdd 
is the strength of the dipolar interactions, 9 is the angle 
between the vector connecting two dipolar particles and 



the dipole orientation, g 
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length, Vt is the trap potential 
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where 7 = — is the trap aspect ratio. It is well known 

that the sign of the dipolar mean-field energy can be 
controUed via the trap aspect ratio. For dipolar conden- 
sates, it is useful to introduce a dimensionless parameter 
that measures the relative strength of the dipolar and s- 

wave interactions Sdd — — — • Chromium atoms posses 

an anomalously large magnetic dipole moment, hence 
e — 0.16. It can be enhanced via Feshbach tuning of the 
scattering length. In Thomas- Fermi regime, dipolar BEC 
is stable as long as —0.5 < Sdd < 1- 

The equation ([ij is an integro-differential equation since 
it has both integrals and derivatives of an unknown wave 
function. Let us focus on the dipolar potential term, 
^dd- There are two main methods to calculate ^dd- 
The first one is based on a Gaussian trial function, ^dd 
can be calculated by employing Fourier transform from 
(r, z)-space to {k±, kz)-space. The Fourier transform of 
the dipole-dipole interaction term in the small cut-off 

distance is given by J"{ Yi-iT'i'^ } = ^ (Scos^a - 1) 
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[Tol . Hence, the convolution theorem leads 

'J'.. = ^-F-i{(3cos2a-l)J-{|*P}} (3) 

where J- and J-"~^ denote the Fourier and the inverse 
Fourier transforms, respectively and a is the angle be- 
tween the momentum k and the dipole direction. Since 
we have assumed that the dipoles are polarized along 

the z-axis, cos^ a = — — ^^2"- The corresponding dipole- 

dipole energy can be calculated using eq. ([3]) and the 
Parseval's theorem 

Em=^ y(3cos2a-l)|^{|vI/|2}|2d3j^ (4) 

where d'^k = 27r k±dk± dk^- The last relation is of im- 
portance in the variational method, which is an efficient 
method to study the stability properties and collective 
excitations of dipolar gas j9l-ll9l|. 

The second method to calculate ^dd was presented by 
Eberlein, Giovanazzi and O'Dell in TF regime ^20|. They 
showed that a parabolic density profile remains an exact 
solution for an harmonically trapped vortex-free dipolar 
condensate in the TF limit. In the TF regime, the mean- 
field dipolar potential, $dc;(r), is given by 

(5) 

where the function /(k) for oblate case, (k > 1), is given 

2 + K^(4-6 '^■^'^ta-ly/;?^ N 

««' 2(1 -.'f^ - 

the condensate aspect ratio. In the absence of dipolar 
interaction, the condensate aspect ratio and the trap as- 
pect ratio match. 

III. A SINGLE VORTEX 

The density profile of a dipolar condensate with a 
straight central vortex line in the TF regime reads (5[ 

where n(r) — when the right hand side is negative and 
/?, -R, and L are variational parameters that describe the 
size of the vortex core, the radial and the axial sizes, re- 
spectively. O'Dell and Eberlein derived the correspond- 
ing energy expression elegantly and performed numeri- 
cal calculations for a prolate trap with 7 = 0.2 and an 
oblate trap with 7 = 5 [H]. They found that the criti- 
cal angular velocity, ilc, above which a vortex state is 
energetically favorable, is increased for the former one 
and decreased for the latter one. Based on the results 
for 7 = 0.2 and 7 = 5, the authors concluded that 17^ is 
decreased due to the dipole-dipole interaction in oblate 



traps, 7 > 1, and increased in prolate traps, 7 < 1. We 
find that this generalization does not hold. The crossover 
between a reduced and an increased critical velocity oc- 
curs at 7 = 7c > 1 rather than 7 = 1 in TF regime. In 
other words, Qc are the same for dipolar and non-dipolar 
condensates at jc- The critical value jc changes slightly 
with Cdd- For ^^Cr atoms with £dd = 0.16, 7c~2.8. For 
Edd = 0.6, 7c « 2.9. Fig-[T] plots fie in units of uj± as a 
function of trap aspect ratio 7 for three different values 
of Edd {£dd = 0, Edd = 0.2, and Sdd = 0.6). As can be seen 
from the figure, f2c increases with Sdd for 7 < 7c, while 
flc decreases with Edd for 7 > 7c. 

So far, we have focused on a central vortex. In [8], based 
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FIG. 1: The critical angular velocity of a condensate for 
£dd = 0, Sdd = 0.2, and Edd ~ 0.6 as a function of trap as- 
pect ratio 7. Qc is measured in units of iu±. Jlc increases 
(decreases) with Sdd for 7 < 7c (7 > 7c), where 7c ~ 2.8 

on two examples, it was concluded that the effect of the 
dipole-dipole interaction is the lowered (raised) preces- 
sion velocity of an off-center straight vortex line in an 
oblate (a prolate) trap. We conclude that the effect of 
the dipole-dipole interaction is the lowered (raised) pre- 
cession velocity of an off-center straight vortex line for 
7 < 7c (7 > 7c)- 

Having found a critical trap aspect ratio in TF limit, let 
us search for the existence of a critical trap aspect ratio 
in Gaussian approximation. 



IV. STABILITY 

In this section, we will explore another critical value 
missed in earlier papers investigating stability properties 
of harmonically trapped dipolar condensates. Consider 
a pure dipolar condensate polarized in the z-direction. 
Note that due to the presence of a Feshbach resonance, a 
practically pure dipolar system can be produced by tun- 
ing the s-wave scattering length to zero. 
A variational method can be employed to obtain an es- 
timate for the instability conditions. We will use the 
notation used in Q. 

1 1 
7^|2 ' • (7) 
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We employ the following trial function 



N 



exp 



where d± — \/h/muj± is the oscillator strength and the 
variational parameters a measures the change in radial 
dimension and k is the condensate aspect ratio. 
The total energy = + E°,, where E°^, are 
dipolar and non-dipolar parts, is given by 
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where E±^ — -N'huj±, the function E{k) 
is a continuous function (k > 0) and 

TO Cdd 
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We investigate the existence of minimum of the total en- 
ergy as a function of the two variational parameters, a 

dE 

and K. The minimum conditions are given by — — = 

OCT 

OE 

and — — = 0. If we substitute the total energy, E'^ , into 

OK 

these equations, we obtain 



TV* 



4 (K-*-i4)-*Ke 



18 



((7 



)F + 2 



11 



^ ' K'^ K* K'^ 



10 



(12) 

The first one enables us to determine the condensate as- 
pect ratio, K, for a given number of particles and trap 
aspect ratio. The second equation gives the information 
how the radial and axial lengths change with k and I. 
To find the critical number of particles, we should also 

calculate the Hessian, A = „ „ — „ ^ 

The local minimum of energy disappears and instabil- 
ity of the condensate occurs when the Hessian is equal to 
zero. The Hessian is too long to write here. Below, we 
will predict the existence of a critical trap aspect ratio. 



Only the first critical point in Gaussian limit was pre- 
dicted in the literature Q. The critical trap aspect ra- 
tios play important roles in understanding of the gen- 
eral structure of dipolar condensates. Firstly, the dipole- 
dipole interaction is positive when I < I* {I > I**) for 
Cdd > {Cdd < 0). Hence, the dipolar condensate will 
be stable at any number of particles as long as Z < Z* for 
Cdd > and I > I** for Cdd < 0. The dipole-dipole in- 
teraction is negative when I > I** {I < I*) for Cdd > 
{Cdd < 0). So, the stability requires TV* < iV* where 
N* is the critical a critical value. The critical number 
A^*, contrary to the case of contact interaction, depends 
strongly on the trap geometry. Let us suppose k — Kc and 
a — Oc when N* — N*. It is interesting to observe that 
N*KiQ.2, Kc~2.6, and tTc~0.68 remain almost the same 
when Z > 3 for Cdd > 0. Fig-Hand Fig-Oplot |iV*| (solid 
curve) and Kc (dashed curve) as a function of trap aspect 
ratio for Cdd > and Cdd < 0, respectively. Amazingly, 
N* and Kc change sharply at I — 0.8 and I ~ 0.7, respec- 
tively for Cdd > and at Z = 2.3 and I = 2.1, respectively 
for Cdd < 0. As can be seen from the figures, there ex- 
ists a minimum for critical number, iV* = 4.64 at I = 0.7 
for Cdd > 0, while no such minimum occurs for Cdd < 0. 
We note that in the region I* < I < I** , the sign of the 
dipole-dipole interaction depends on Cdd and I. 
Secondly, the condensate aspect ratio is equal to the 
trap aspect ratio, k — I, only if I = I* or Z = I**. Only 
in the region I* < I < I** , k is bigger (smaller) than 
I for Cdd > {Cdd < 0). As N* is increased, k in- 
creases (decreases) when I* < I < I** {I < I* or: I > I**) 
for Cdd > 0. Conversely, as iV* is increased, k de- 
creases (increases) when I* < I < I** {I < I* oy I > I**) 
for Cdd < 0. The dipole-dipole interaction stretches 
(squeezes) the cloud both in the radial and the axial 
directions compared to the non-dipolar condensate if 
I < I* for Cdd > {Cdd < 0) and if I > I** for Cdd < 
{Cdd > 0). 

It is remarkable to note that the existence of two crit- 
ical trap aspect ratios, 0.44 and 2.30 (in our nota- 
tion), for a metastable TF solution outside of the regime 
—0.5 < Edd < 1 was discussed [22,]. It is interesting to ob- 
serve that these critical values in TF approximation are 
very close I* and I** in Gaussian approximation (fT3|) . It 
was shown that collapse via scaling oscillations is sup- 
pressed in TF limit if Z > 2.30 for Sdd — > — 00 and if 
Z < 0.44 for Edd 00 [2^ . 

In the following two subsections, we will show that the 
critical trap aspect ratios are changed by an additional 
anharmonic potential 
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Z64 



(14) 



A. Critical Points 

There exists two critical values for the trap aspect ratio 
Z* = 0.43 , I** = 2.52 . (13) 



where A_l and Xz are small dimensionless parameters 
characterizing the strength of the quartic potential in ra- 
dial and axial directions, respectively. This potential has 
been achieved experimentally by superimposing a blue 
detuned laser beam with the Gaussian profile to the mag- 
netic trap [2l| . 
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FIG. 2: The critical number of atoms, A'^^, (solid) and the 
critical value of the condensate aspect ratio, Kc, (dotted) as 
a function of the trap aspect ratio, I, for a purely dipolar 
condensate with Cdd > 0. 




FIG. 3: The critical number of atoms, |A^*|, (solid) and the 
critical value of the condensate aspect ratio, Kc, (dotted) as 
a function of the trap aspect ratio, /, for a purely dipolar 
condensate with Cdd < 0. 



B. Axially Quartic Potential 

We will now investigate the critical trap aspect ra- 
tios in the presence of anharmonic potential. Suppose 
first that anharmonic potential is in the axial direction, 
Aj_ = 0. The selection of the proper form of trial func- 
tions is very important in the variational approach. We 
will choose a trial function which reduces to the ground 
state solution of the Schrodinger equation in the nonin- 
teracting limit. Hence, we employ the trial function up 
to the first order of X, 



4" = ■00 exp 
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I iV 9 
where ipo = \ 3,^ 3 3 — (1 + — A^) is the normahza- 

tion constant. If we evaluate the energy expression, 
-£^nd + ^dd- weget 
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where E± = —NhiL!±, E^^i and are the dipole-dipole 

and non-dipole energies for purely harmonic trap ^ and 
G^(/c) = 14 - SS/t^ + 52k^ + 7k*'. The minimums of the 
variational parameters a and k can be calculated from 
these energy expressions. Then critical aspect ratios can 
be found. We find that they are decreased with A^. For 
example. 



?L=o = 0-43; ^L=o.05 = 0.41; 



;j^^o.i = 0.39. (17) 



11%^ = 2.52, 11%^,^^ = 2.43 lll^^,, = 2.34 . (18) 

The critical number of particles, doesn't change ap- 
preciably with Xz- When I > I**, the critical values of 
the condensate aspect ratio, Kc, is decreased with Xz for 
Cdd > 0. For example, if Z = 10 then Kc,a^=o = 2.60, 
Kc,A^=o.05 = 2.49 and Kc,\^=o.i — 2.40. However, CTc ~ 
0.68 slightly increases with Xz- 



C. Radially Quartic Potential 

Suppose now that anharmonic potential is in the radial 
direction, A^ = 0. We will choose our trial function up 
to the fist order of Aj^ 



4* = ■00 exp 
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(r' + ^ + X^{r^ + 



4^24 



)) 



(19) 

where ipo — \ / , ^ — 5~(1 H — ^^-l) is the normalization 

constant. The corresponding energy, E'^^ = E^^^ + E^j^ , 
is given by 
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where G_l(k) = 112 -h 45^2 - 258k^ + 56k^ 
From the energy expression, we can find the critical trap 
aspect ratios. The first critical point changes very slowly 
with A_L in contrast to the case of axially anharmonic 
potential (fT7|) . 
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The second critical point is increased with Aj^. For ex- 
ample, 



7* 



2.52, I 



1 'Ai=0.05 



2.56 l\ 



2.61 



(22) 



Nc and Cc doesn't change appreciably with Aj^. When 
I > I**, Kc is increased due to the quartic potential. As an 
example, if I = 10, then Kc.x^=o ~ 2.6, Kc,x±=a.Q5 — 2.7 
and Kc,A , =0.1 = 2.8. 
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In this paper, we have shown that there exists some crit- 
ical trap aspect ratios for a dipolar condensate determin- 
ing the general structure of the condensate. Classifying 
the trap as prolate and oblate is not helpful. It was gener- 
ally believed that critical angular velocity, ilc, decreases 
(increases) due to dipolar interaction in an oblate (a pro- 
late) trap. Here, we have discussed that there exists a 
critical trap aspect ratio 7c~2.8 above (below) which flc 
is increased (decreased) due to dipolar interaction. We 



have also shown that the same is true for the preces- 
sion velocity of an off-axis vortex. Furthermore, we have 
proved that there are two critical trap aspect ratios for 
a purely dipolar gas in Gaussian limit. The stability di- 
agram depends crucially on the critical trap ratios. We 
have shown that the critical points can be changed in the 
presence of quartic potential. 

Helpful discussions with Z. Oztas are gratefully acknowl- 
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